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Abstract 

We study the large-order asymptotic behavior of the perturbation series for 
short-distance coefficients in the NRQCD factorization formulas for the de- 
cays J/tp — > e + e~ and r\ c — > 77. The short-distance coefficients of the leading 
matrix elements are calculated to all orders in the large- Nf limit. We find that 
there is a universal Borel resummable renormalon associated with the cancel- 
lation of the Coulomb singularity in the short-distance coefficients. We verify 
that the ambiguities in the short-distance coefficients from the first infrared 
renormalon are canceled by ambiguities in the nonperturbative NRQCD ma- 
trix elements that contribute through relative order v 2 . Our results are used 
to estimate the coefficients of higher order radiative corrections in the decay 
rates for J/tp — ► e + e~, r\ c — ► 77, and J/tp — ► 777. 
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Typeset using REVIgX 



I. INTRODUCTION 



The annihilation decay of heavy quarkonium involves several distance scales. The annihi- 
lation of the heavy-quark pair is a short-distance process occurring at length scales of order 
1/m, where m is the mass of the heavy quark. The probability that the heavy-quark pair will 
be close enough to annihilate is determined by the wavefunction of the heavy quarkonium, 
which is dominated by distances of order 1/mv or larger, where v is the typical relative 
velocity of the heavy quark. The nonrelativistic QCD (NRQCD) factorization formalism 
HIIJ] is a systematic framework for separating the effects of the scale m from effects of lower 
momentum scales. The annihilation decay rate is factored into short-distance coefficients 
and long-distance NRQCD matrix elements. The short- distance coefficients can be calcu- 
lated using perturbative QCD as a power series in the strong coupling constant a s {m) at 
the scale of the heavy quark mass. The matrix elements are nonperturbative, but scale in 
a definite way with v. The NRQCD factorization formalism therefore organizes the decay 
rate into a double expansion in powers of a s and v. 

In principle, the short- distance coefficient for each matrix element can be calculated 
to any desired order in a s by matching the NRQCD effective theory with full QCD. One 
might hope that any desired precision could be obtained by carrying out the matching 
calculations to sufficiently high order in a s . However, large numbers may arise in high order 
calculations through the factorial growth of the coefficients in the perturbation series. A 
factorially growing series, when truncated at a high order, may produce a large error even if 
the coupling constant is small. This may limit the possibility of obtaining sufficiently precise 
results by performing the perturbative calculations to higher orders. 

One can investigate the asymptotic behavior of a series f(a s ) by studying its Borel trans- 
form f(t). A singularity in the Borel parameter t is called a renormalon ||. A renormalon 
gives rise to n\ growth of the coefficients in the power series for f(a s ) so that its radius of 
convergence is 0. If the renormalon is on the negative real axis, Borel resummation can be 
used to recover the function uniquely from the divergent series. If the renormalon is on the 
positive real axis, the sum of the divergent series is not unique and there is an ambiguity in 
the function f(a s ). 

The renormalons that are responsible for the factorial growth of coefficients in the per- 
turbation series for QCD have been extensively studied. Important applications include 
the operator product expansion for e + e~ annihilation and r decay [f§-0] and heavy quark 
effective theory The large- Nf limit of QCD has been developed by Beneke || into 



a useful tool for the explicit study of renormalons. In this paper, we apply this tool to the 
NRQCD factorization formulas for S'-wave quarkonium states, focusing on electromagnetic 
annihilation decays. 

The rest of this paper is organized as follows. In section II, we calculate the short- 
distance coefficients for the electromagnetic annihilation decays J /if) — > e + e~ and r] c — > 
77 at tree level. We present a new method for calculating the tree-level short-distance 
coefficients which involves expanding nonlocal QCD matrix elements in terms of local QCD 
matrix elements, and then using a Foldy-Wouthuysen-Tani transformation to express them 
in terms of NRQCD matrix elements. In section III, we calculate the Borel transforms 
of the short- distance coefficients of the leading matrix elements in the large- Nf limit. In 
section IV, we study the renormalons in the short- distance coefficients. We show that there 
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is a universal Borel resummable renormalon related to the cancellation of the Coulomb 
singularity. We demonstrate that the ambiguities from the first infrared renormalon are 
canceled by ambiguities in the matrix elements that contribute through relative order v 2 . 
We also deduce the leading renormalons for the decay J /if) — > 777. Finally, in section V, 
we use our calculations to estimate higher order radiative corrections to the decay rates for 
J /if) — > e + e~ , T] c — > 77, and J /if) —* 777. 



II. SHORT DISTANCE COEFFICIENTS AT TREE LEVEL 

In this section, we review the NRQCD factorization formalism as it applies to electro- 
magnetic annihilation decays. We then calculate tree-level short-distance coefficients for the 
decays J/if> — > e + e _ and rj c — > 77 and for several quarkonium-to- vacuum matrix elements. 



A. NRQCD Factorization Formalism 

In the NRQCD factorization formalism, the annihilation decay rate for the charmonium 
state H is written in the factorized form [0J 

r W = oTF~ E C mn(H\O mn \H) , (1) 



2M H 



mn 



where Mh is the mass of the state H. The matrix elements (H\O mn \H) are expectation 
values in the quarkonium state H of local 4-fermion operators that have the structure 

O mn = ^JC mX , (2) 

where if) and \ are the field operators for the heavy quark and antiquark in NRQCD, K n 
and K m are products of a color matrix (1 or T a ), a spin matrix (1 or a 1 ), and a polynomial 
in the gauge-covariant derivative D and in the field strengths E and B. It is convenient 
to define the states \H) = \H(P = 0)) in (fl|) so that they have the standard relativistic 
normalization: 

'H{P') H(P)) = 2E P (2tt) 3 5 3 (P - P') , (3) 



where Ep = yM| + P 2 . The coefficients C mn take into account the effects of short distances 
of order 1/m, and they can therefore be calculated as perturbation series in the QCD coupling 
constant a s (m). 

For electromagnetic annihilation decays, the factorization formula can be simplified. 
Because there are no hadrons in the final state, one can insert a projection |0)(0| onto the 
vacuum state between \ and i n the operator (|J). The factorization formula ([I]) then 
becomes 

T(H - EM) = -L Y,Cmn{H\if)*JC mX \0){0\Am) , (4) 

lM H mn 

where EM represents an electromagnetic final state consisting of photons or e + e~ pairs. The 
matrix elements vanish unless the operators are color singlets, so K n and K m are products 
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of a spin matrix (1 or a 1 ) and a polynomial in D, E, and B. The short-distance coefficients 
C mn in (|4]) include the effects of integrating over the phase space of the electrons, positrons 
and photons in the final state. The NRQCD factorization can actually be carried out before 
the phase space integration at the level of the T-matrix element: 

Th^m = Y, C m(®\x ] K m ^\H) • (5) 

m 

The velocity-scaling rules of NRQCD can be used to determine how the matrix ele- 
ments in (H) scale with the typical relative velocity v of the heavy quark. For decays of 
the S-wave states J /if) and rj c , the leading contributions are from NRQCD matrix elements 
{Q\x^ cnp\J / if)) an d (O|x^0|?7c}j which scale like v 3 ^ 2 . The next most important matrix el- 
ements are (0|;\^D 2 cn/>| J/ip) and (0|x^D 2 ^|?7c), which are suppressed by v 2 and therefore 
represent relativistic corrections. The matrix elements are related by the approximate heavy- 
quark spin symmetry of NRQCD: 

(OlxWl JMe)> = e(0\ X ^\Vc) (l + 0{v 2 )) , (6) 
(Olx+DV^JMe)) = 6(01x^1^) (l + 0(v 2 )) , (7) 

where e is the polarization vector of the J /if} which satisfies e 2 = 1. Thus there are only 3 
independent matrix elements that contribute to S'-wave decays through relative order v 2 . 

The short-distance coefficients in the factorization formula (||) can be determined by 
the covariant projection method \T2\ or by the threshold expansion method flT3| |. In the 



threshold expansion method, the coefficients are determined by matching perturbative T- 
matrix elements for cc annihilation. Let |cc) = |c(q, £)c(— q, rf)) represent a state that 
consists of a c and a c with spatial momenta ±q in the charmonium rest frame and with 
spin and color states specified by spinors £ and r\ that are normalized to £t£ = rfr) = 1. The 
standard relativistic normalization is 

c(qi,Oe(^.V)|c(qi,0c(qa,»7)) = ±E qi E q2 (27r) 6 5 3 ( qi - q'^ 3 ^ - q' 2 ) , (8) 

where E„ = y/m 2 + q 2 . In the threshold expansion method, the matching condition is 



■>EM 



pQCD 



Eft.Mx'K.*!^ • (9) 



where 7^ £ ^em is the T-matrix element for the annihilation of the cc into the electromagnetic 
final state EM. To carry out the matching procedure, the left side of @ is calculated using 
perturbation theory in full QCD, and then expanded in powers of q. The matrix elements on 
the right side of are calculated using perturbation theory in NRQCD, and then expanded 
in powers of q. The short-distance coefficients C m are obtained by matching the terms in 
the expansions in q order by order in a s . To determine all the coefficients C m , it is also 
necessary to consider the annihilation of cc states with general momenta qi and q2 and the 
annihilation of higher Fock states like |ccg). At tree level, some of these coefficients can be 
determined more easily using an alternative method that involves a Foldy-Wouthuysen-Tani 
transformation. This method will be presented in the next subsections. 



4 



B. Leptonic Decay of J/ip 



The decay J/ip — > e + e~ proceeds through the annihilation of the J/ip into a virtual 
photon. The Feynman diagram for cc — > e + e~ is shown in Fig. 1. The T-matrix element for 
the decay of J/ip into an electron and a positron with momenta k± and k 2 can be written as 

Tj,^ e+e - = -9^u(h)^v(k 2 ) (0\*-f*\ JM , (10) 

where Q = +2/3 is the electric charge of the charm quark and is the mass of the J/ip. 
The factor of l/M^ comes from the propagator of the virtual photon. The matrix element 
in (|1(]) is that of a local QCD operator between the J/ip and the vacuum state. Here and 
below, the argument of a local operator is assumed to be x = if it is not given explicitly. 
We obtain the decay width by squaring the T-matrix element, summing over the spins of 
the e + and e~ , and integrating over the phase space of the e + e~ pair: 

r(J/^ - e+e-) = -^L_ {-g, u ) (J/V^t^O) (0|* 7 ^l W> • (11) 

The QCD matrix element (0| 1 4 r 'y A *^E r | J/if>) in (pi) includes effects of both short distances 
of order 1/mor smaller and long distances of order 1/mv or larger. We can separate these 
effects by using the NRQCD factorization formalism to express the QCD matrix element as 
a sum of NRQCD matrix elements multiplied by short-distance coefficients. At tree level, 



this is most easily accomplished by a Foldy-Wouthuysen-Tani (FWT) transformation flU 
This transformation can be constructed order by order in v so that off-diagonal terms in 
the Hamiltonian that couple the heavy quark and antiquark are suppressed to any desired 
order in v. Up to corrections of order i> 3 , the FWT transformation is 



¥(*) = exp -i-y-D I" , (12) 




where D = V — ig s A is the covariant derivative. Inserting QT2D into the matrix element 
(0| 1 l , 7 /i ^/| J/fp), expanding to second order in D/m, and keeping only the annihilation terms 
that involve an d ip, we obtain 



(0|*7 M #|J/V>) = -g^ 



(o| x VV|JM + ^(oix^AA' + ^-A)^!^) 



(13) 



Up to this point, the only property of the J/ip that we have used is that it has the quantum 
numbers J PC = 1 that allow it to annihilate into a virtual photon. According to the 
velocity-scaling rules of NRQCD, the first term on the right side of ([13]) scales like v 3 / 2 . 
The last term in ( |T3| ) is a symmetric tensor in the indices of the two covariant derivatives. 
The trace of this tensor contributes at relative order v 2 . The traceless part receives a 
contribution from the D-wave component of the |cc) Fock state of the J/ip and is suppressed 
by v 4 . Keeping only the terms in (|l^) that contribute through relative order v 2 , this matrix 
element reads 

(o|*y*|J/v> = -<r (0|xVv| J/4,) + -^(OI^DVVIJM 



(14) 
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Inserting ( ^4|) into (|TT]) , the decay width can be written as 



4nQ 2 a 2 
3M| 



e^e 



(olx^viJM + ^(ol^DV^lJ/V) 



(15) 



This is the NRQCD factorization formula for the decay width, including the first relativistic 
corrections and with the short- distance coefficients calculated to tree level. 

In Ref. 0, the factorization formula for the decay rate of J/ip — > e + e~ was derived by 
matching cc scattering amplitudes. The T-matrix element for cc — > cc through an interme- 
diate e + e~ state was matched with the T-matrix element for cc — > cc through operators of 
the form ^ K mX \$) <0|x j 'K n ip in NRQCD. Using our relativistic normalization of states and 
keeping only terms through relative order v 2 , the result is 



r(J/V 



nQ 2 a 2 
3M^m 2 



e + e 



(0\ X *<ri/>\J/1>) + 



2 (ol^DV^I J/V) 12 



3m 2 



(16) 



As pointed out by Maksymyk |jl5| , the consistency of (fl5|) and (pl|) through relative order 
v 2 requires the identity 



(0|x f DV^|J/^) = -m(M i ,-2m)(0\x j (Tip\J/*P} (l + 0{v 2 )) 



(17) 



This identity was recently derived by Gremm and Kapustin [ IB] , and follows simply from 
the equations of motion of NRQCD. It allows factors of m in a factorization formula to be 
traded for factors of M^: 



2m 



1 



1 



2m 2 e- (0|xWIW(c)) 



+ 0(v 4 



(18) 



The identity ( |18[) can also be used to relate our NRQCD matrix elements, which are de- 
fined using the standard relativistic normalizations for the quarkonium states, to the matrix 
elements of Ref. [II], which were defined using the standard nonrelativistic normalizations: 



(0| X f /C>l^) = v^M^OIxXVWbbl- 
Using (|T8| ) and (|19"D , we obtain the relation 



(19) 



/BBL 



Am 2 



(OI^DVWMbbl 



0(v 4 



(20) 



Inserting (|20| ) into our factorization formula ( [T5| ) and using ([T8|) to eliminate in favor of 
m, we recover the result given in Ref. Ill] up to corrections that are higher order in v : 



2nQ 2 a 2 
3m 2 



e e 



(oix^viJM 



BBL 



2 (OI^DV^J/^bbl' 2 



3m 



(21) 
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C. Decay of r/ c into Photons 



If a quarkonium state H is even under charge conjugation, it can decay into two photons 
at second order in the electromagnetic interaction. The Feynman diagrams for cc — > 77 are 
shown in Fig. 2. If the photons have momenta k\ and k 2 and polarization vectors e\ and e 2 , 
the T-matrix element for the decay is 

Tu^n = ~16Q 2 e 2 e ltl e 2u J d 4 z e l ^- k2) z 

x (0| [ f {z) Y S (z, -z) 7 "tf {-z) + * (-z) Y S {-z, z) (z) ] \H) . (22) 

The factor of 16 arises from taking the integration variable z to be half the separation of the 
interaction points. This expression involves the full propagator S(x,y) of the heavy quark, 
which satisfies the equation 



(ifl- m) x S(x, y) = 5 4 (x - y) , 



(23) 



where D^ = + ig s A^. 

The phase factor in ( p2]) and the behavior of the quark propagator combine to give the 
integral over z support only in the region where the components of z M are of order of 1/m. 
Since this is much smaller than the size 1/mv of the quarkonium, the operator in ( [22] ) can 
be expanded in powers of z^. It is convenient to carry out this expansion using the radial 
gauge 0: 



z»AJz) = . 



(24) 



By Taylor expanding this gauge condition around z — 0, one can show that derivatives of 
at the origin can be expressed in terms of gauge-invariant operators. In particular, we have 
Afj,(0) = and <9 M A,(0) = F llu (0)/2. The quark field in the radial gauge can be expanded as 



*(*) 



tf(0) + z»D^(0) + -z' l z u D lt D v ^(0) + 



(25) 



By solving (|23f) , we find that the quark propagator can be expanded as 
S (z, —z) = Sp(2z) + J d 4 x Sp (z — x) g s 4-( x ) Sp (x + z) + • • • 



-2ip-z 



1 



$ — m 



(p 2 — m 2 ) 2 



Y(jt-m)Y + 



(26) 



where Sp(x) — L e tp ' x /(p'—m) is the Feynman propagator for a free quark. We have 
introduced an abbreviated notation for the integral over a four-momentum: 



(27) 



Inserting (|25|) and ( pB"l) into (|22]) , expanding up to operators of dimension 5, and integrating 
over z, we obtain 
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Th-tt, = ~Q 2 e 2 e lfl e 2v (0|* -(r A T^ + 2m^ 



■ r 2 — m 2 

+ 4 7T~ -I^—i ( r A rt Al/ + mlT) A, 
ar a r z — m z v y 

2 or a or p r z — m z v ' 
~ \ p _ 1 m2)2 ( r> - mT^ ) ig s F a , } 9\H) . (28) 

We can set r = (k± — k<i)j2 and r 2 = —Mfj/A after taking the partial derivatives. We have 
used the notation 

pMiM2---Mjv _ ^fii . . . ^Miv _|_ ^mjv . . . y2 _ ^29) 

To further simplify the T-matrix element, it is convenient to choose the photon polarization 
vectors to be space-like in the rest frame of the H. In that frame, we have r = 0, r • e% = 0, 
and r • e 2 = 0. The FWT transformation ( |T2"D can be used to reduce T to a sum of NRQCD 
matrix elements. Keeping all terms through second order in D/m, we obtain 

T »- = - mB^ ^ { ^ + D "" + D >"' + tSt£s d ™°" 

(Mfj + 12m 2 ) e ijk r k 2 64e^' fc r fc r"V n 
+ 2m 2 (M^ + 4m 2 ) D ~ [M 2 H + 4m 2 ) 2 m n 



(M% - Am 2 ) e ijk r k 
h 2m 2 (M 2 H + Am 2 ) " ' °" 
4m 



M|f + 4m 2 



g s {E l a 3 + £V ! - 5 y E • &) } i/>\H) . (30) 



To determine the coefficients of the matrix elements involving B and E using the threshold 
expansion method, we would have to match scattering amplitudes for ccg — > 77. 

Up to this point, the only property of the quarkonium state H that we have used is that 
its charge conjugation quantum number is C = +1, so that it can decay into two photons. 
Thus the formula ( |30"D applies equally well to the r) c , the Xco, and the Xc2- We now specialize 
to the 77c, which has quantum numbers J PC = h . The parity conservation of NRQCD 
implies that the matrix elements of operators linear in D or linear in E are zero. Using the 
rotational invariance of NRQCD, the remaining matrix elements in (|3~0"D can be reduced to 
scalars: 

_ . AQ 2 e 2 M, c . 
7 " M 2 +4m 2 eiX62 - r 

f / > + , . x 3 K + 16m 2 M 2 + 144m 4 , . . _ , ( x 

x {<°i^i*> + \«mi f^f <0|x D * w 

M 2 - Am 2 , 1 

+ 2 m »(< + 4m») ("lAB.^fe)). (31) 
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According to the velocity-scaling rules, the matrix element (0|x^|?7c) scales like v 3 ! 2 . The 
matrix element (0|x^D 2 ^|^ c ) represents a relativistic correction that is suppressed by a 
power of v 2 . The matrix element (0|x^sB • <Tip\r] c ) is suppressed by v 7 ! 2 . This suppression 
factor can be deduced using the methods described in Ref. 0. In Coulomb gauge, the 
leading contribution to this matrix element comes from a \ccg) Fock state, with the cc in 
a color-octet 3 Si state. This Fock state is dominated by dynamical gluons with momenta 
of order mv and has a probability of order v 3 . Relative to (0|x^l^c), the matrix element 
(0|x^B • (Tijj\r] c ) scales like (mv) 2 v 3 ^ 2 , with the factor of v 3 ^ 2 coming from the amplitude for 
the ccg Fock state and the factor of (mv) 2 coming from the dimension of the operator B. 
Thus the suppression factor is v 7 ! 2 . Keeping only the terms through relative order v 2 , the 
T-matrix element reduces to 



T 



—I 



4Q 2 e 2 M rjc 
M 2 + 4m 2 



ei x e 2 ■ r 



f , , t , v 3M « + 16m 2 M n 2 + 144m 4 , 
X { <°l^> + W^ 2 ) 2 ^ ^ 



(32) 



We obtain the decay rate by squaring the T-matrix element, summing over the polariza- 
tions of the photons, and integrating over their phase space: 



r(?7c -> 77) 



IQttQW M Vc 
(M 2 c + 4m 2 ) 2 



x 



. + , . . 3M 4 + 16m 2 M 2 + 144m 4 , , , „ , , x 1 



(33) 



The factorization formula ( |33"D can be greatly simplified by using the identity analogous to 
(H) for the r} c : 



2m 



1 (Olx+DVl^) 



2m 2 (0\ X ^\Vc) 



+ 0(v 4 ) . 



(34) 



Eliminating the quark mass from the coefficient of the leading matrix element in (|33|) and 
keeping only terms through relative order v 2 , we get 



r(?7c -> 77) 



4vrg 4 a 2 



(o\xty\vc) + 



6m 2 



IX 



(35) 



The matrix elements can be expressed in terms of those introduced in Ref. [TJ by using ( |T9"D . 
Using (|34] ) to eliminate the remaining factors of M r?c in favor of the quark mass, we recover 
the factorization formula derived in Ref. |j| up to corrections that are higher order in v: 



F(Vc -> 77) 



27rg 4 a 2 



m z 



(olxVl^c) 



BBL 



+ 



3m 2 



(Ol^D^I^BBL 



(36) 
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D. Quarkonium-to- Vacuum Matrix Elements 



In order to deduce the renormalon ambiguities in the NRQCD matrix elements in the 
factorization formulas fll5 ) and fl35|), it will prove useful to also consider quarkonium-to- 
vacuum matrix elements for all the QCD operators with dimension 3. A complete set 
of these operators is the scalar density tyty, the vector current 1 J r 7 /i \E r , the tensor current 
tya^ty, the axial vector current ^7^75^/, and the pseudoscalar density \j/7 5 \I/. Using the 
FWT transformation ([12]), we expand the matrix elements of each of these operators in 
terms of NRQCD matrix elements. Keeping terms up to second order in D, the matrix 
elements for a general quarkonium state H are 



(OlMif) = — (OMD • <Tip\H) 
m 



(0|^ 7 ^|#) = 8* <0|x + 



1 



1 + 



4m 2 
1 ( 

2m? 



DjDi) a J 



D 2 + g s B a 



+ 5^ — <0|x f (D x tryi/>\H) 



(0\^a^\H) = 8^5 uj — e ijk (0\ X ^D k iP\H} - z(<^° ST - 8 u0 5*) 

rf' (2DV {D i D j • l)jl),)o J +2g s B i 



x (0|x 



Am 2 



(0\^-f 5 ^\H) 



(0\ X ty\H) . 



(37) 
(38) 

(39) 

#ff>, (40) 
(41) 



We now specialize to the J/ip and the r] c . For the J/i/j, the only operators that have 
matrix elements of order t> 3 / 2 are the vector current, whose matrix element is given in (|I"4"D, 
and the tensor current. Keeping only term through relative order v 2 , the matrix element of 
the tensor current is 



Itf^tflJ/VO = -i{8^ 8 vi - 8 vU 8^ 



(0| X VV|J/V) + ^(olx+DVVIJ/v) 

6m A 



(42) 



For the r/ c , the operators that have matrix elements of order v 3 ^ 2 are the axial vector current 
and the pseudoscalar density. Keeping only term through relative order v 2 , these matrix 
elements are 



(0|*7" T5 *|77 C ) = 8"° 



(o|xVk) + 



2m [ 



:{Q\ X ^\Vc) 



(0\^\ Vc ) = -(OlxVl^c) 



(43) 
(44) 



III. SHORT-DISTANCE COEFFICIENTS IN THE LARGE-iV F LIMIT 



In this section, we briefly review the Borel resummation method and the standard pre- 
scription for calculating the Borel transform of a perturbation series in QCD in the large- Nf 
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limit. We apply this prescription to the pole mass to illustrate the renormalization of 
the Borel transform. We then use the threshold expansion method to calculate the Borel 
transforms of the short- distance coefficients of the leading matrix elements in the NRQCD 
factorization formulas for (O^^l J/ip), (0\^a^^\J/ij;), (01*7^75*!^), (0|^7 5 *|r/ c ), and 
the decay amplitude for r\ c — > 77- 



A. Borel Resummation 

If a function f(a s ) has the power series expansion 

00 

fM = a « a * ' ( 45 ) 



n=Q 



its Borel transform is defined by 



00 1 

f(t) = a 8(t) + T^TJj antn ~ l ■ (46) 



n=l 



Obviously, the power series ([46]) for the Borel transform has better convergence properties 
than that of /. The original function / can be recovered from its Borel transform f(t) by 
the inverse Borel transformation: 

f(a s ) = / dt e-^* f(t) . (47) 



If the integral fl47|) is sufficiently convergent, it defines a unique function with the power 
series expansion (|4"5|) and the series is called Borel resummable. 

A singularity in the Borel parameter t is called a renormalon ||. The contribution to 
f(a s ) from a pole of f(t) at t — t* with residue R* can be obtained by replacing f(t) in fl4"T|) 
by R*/(t — t*). Expanding as a power series in t and then integrating, we get 

/(a.) ~ -JTX>-1)! (^)" . (48) 

Thus a renormalon gives rise to n\ growth of the coefficients in the power series for f(a s ). 

If the renormalon is on the negative real axis, this contribution to the series is Borel 
resummable. If the renormalon is on the positive real axis, there is an ambiguity in the 
function f(a s ). The possible prescriptions for the inverse Borel transform include deforming 



the integration contour in ( 47|) into the complex plane so that it runs above the pole or below 
the pole. We also can take a linear combination of these two contours, such as the principle 
value. These prescriptions all differ by amounts that are proportional to the residue of the 
integrand of ( f4"7j ) at the pole. If the pole in f(t) is at the point t* and has residue R*, the 
ambiguity in f(a s ) has the form 

Af(a s ) = K (27r R*) e ~**/ as , (49) 

where K is an arbitrary constant. If a s is the running coupling constant of QCD, its scale 
dependence is approximately 
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where (3q = (33 — 2iV/)/127r. Inserting this into (49), we see that the renormalon ambiguity 
can be approximated by 



Af « K(2wR*) 



^-QCD \ 

. Q 2 ; 



A** 



(51) 



Thus the ambiguity is suppressed by a power of Q 2 - 

To study the renormalons in NRQCD factorization formulas, we need to calculate the 
short- distance coefficients to all orders in a s . An exact calculation of the coefficients to all 
orders would be impossible. However the calculation is tractable in the large-iV/ limit of 
QCD, where Nf is the number of light quarks ||. The diagrams that dominate in this limit 
are obtained by inserting a chain of 1-loop quark bubbles into the gluon propagator in the 
lowest order diagrams that contain a gluon. All other diagrams are suppressed by a factor 
of 1/Nf. QCD is no longer an asymptotically free theory when Nf — ► oo, so it is more 
appropriate to consider the limit Nf —>■ — oo. To be explicit, we define the large- Nf limit to 
be the limit Nf — > — oo, a s — > 0, with Nfa s fixed. 

The insertion of a light quark bubble renormalized by minimal subtraction into the 
propagator of a gluon of momentum I gives a multiplicative factor 



- p a s 



9(e) 



-l 2 -ie 



l 2 + ie 



(52) 



where /3q = —Nf/6ir is the first coefficient of the QCD beta function in the large- Nf limit 

6(47r) e r(l + e) T 2 (2 - e) 



T(4 - 2e) 



(53) 



and 



A = -7 + ln(47r) + - + C 

o 



(54) 



The constant C depends on the renormalization scheme and has the value C = —5/3 in the 
MS scheme. A chain of bubble diagrams forms a geometric series. The Borel transform for 
this series multiplied by a s can be obtained by making the following substitution for the 
gluon propagator in covariant gauge: 



O'.s 



l 2 + ie 



+ (i-0 



exp < u 



9[£) 



exp < u 



9(e) 




,e(d/du' 



(-l 2 -ie) 2 




1 fiu 



-I 2 - ie) l + u ' 



-l 2 - ie) 2 



+ m v terms, 



(55) 
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where u = fat and t is the Borel parameter. The diagrams for electromagnetic decays 
considered in this paper are QED diagrams up to color factors. The contributions from 
the \H V terms in (|5^) cancel by the Ward identities of QED. The calculations can therefore 
be greatly simplified by using a propagator whose Lorentz structure is that of the gluon 
propagator in the Feynman gauge. The substitution ([55]) then reduces to 



P + ie 



Je 



E 

P c 



fj 



111' 



'-l 2 -ieY+ u 



(56) 



where the action of the operator _T t on a function f(u) is defined by 



Jef{u) 



exp < u 



9\ e ) e(d/du') _ ^_ 



dv! 



-(A-C)u' 



/(«0 



(57) 



This operator is defined so that it reduces to the identity operator in the limit e — ^ if it 
acts on a function with no poles in e. 

In calculating the Borel transform of a sum of bubble chain diagrams, the extra factor 
l/(—l 2 ) u introduced in the gluon propagator can introduce divergences in the loop integral 
that do not appear in the one-loop diagram, which corresponds to u — 0. Below certain 
critical values of u, new power ultraviolet divergences appear. At such a value u*, the 
loop integral has a logarithmic divergence, which appears as a pole in the Borel transform 
at u = u*. This pole is an ultraviolet renormalon. Similarly, there are power infrared 
divergences that appear for u above certain critical values u* . At u = u*, the loop integral 
is logarithmically infrared divergent. This logarithmic divergence appears as a pole in the 



Borel transform at u 



u 



and this is an infrared renormalon. 



B. Pole mass and renormalization constant 

We first calculate the Borel transform of the pole mass and the wave function renormal- 
ization constant of the heavy quark. The self-energy for a heavy quark with momentum p can 
be expressed as a function of the coupling constant a s , the matrix ^, and the bare mass m . 
We denote the self-energy by S(^, a s ), suppressing the dependence on m . In perturbation 
theory, the pole mass m po i e (a s ) is defined by the pole in the heavy quark propagator: 



m po i e (a s ; 



m + E(^, a s 



(58) 



rf=m polc (a s ) 



One can solve the equation for m po i e order by order in a s . The on-shell wavefunction 
renormalization constant Z(a s ) is given by Z _1 = 1 — 5Z, where 



SZ{a s ) 



d $ 



(59) 



="Jpol c (a s ) 



A one-loop Feynman diagram for the heavy-quark self-energy with a chain of light quark 
bubbles inserted is shown in Fig. 3. Using the substitution (0), the Borel transform of the 
sum of the diagrams in Fig. 3 is given by 
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where the operator J e is defined in 
dimensional regularization in D = 4 




(60) 



To regularize the ultraviolet divergence, we use 
2e dimensions. The integration measure in ( |60|) is 



tZ^Z 



(27T) 



(61) 



The e-dependent factor in front of the integral has been inserted so that renormalization by 
minimal subtraction will give the same renormalization scheme parametrized by the constant 
C used in (52|). The integral in (|60|) can be reduced to an integral over a Feynman parameter. 
For the Borel transforms of the pole mass in fl5E|) and the wavefunction renormalization 
constant in (|59|), the integral can be expressed in terms of gamma functions: 



m po i e (t) = m 



5Z(t) 



6(t) + 
-2e 



2e 



,(7-5/3)6 



3tt 



Je 



u+e 



e c rriQ 



r(M + e)r(3-2M-2e) 
[l-2u- 2e)T(3 - u - 2e) 



( 7 -5/3)e 



3tt 



e c m, 



u+e 



0, 



1 +M)r(M + e)r(3-2«-2e) 
(1 -2u-2e)r(3-u-2e) 



(62) 
(63) 



The expressions fl62"|) and (|63|) are ultraviolet divergent and require renormalization. A 
renormalization procedure for the large- Nf limit has been derived by Palanques-Mestre and 
Pascual jl8| and by Broadhurst |19|, and applied to the renormalization of the Borel trans- 
form by Beneke and Braun ||. We can simplify their discussion since the renormalization 
of the light quark bubble diagram has already been carried out in fl52|). The ultraviolet 
divergences in (|S2j) and fl63|) appear as poles in e in the coefficients obtained by expanding 
these expressions as power series in u. The poles in e arise from the power series expansion 
of the factor T(u + e). In the minimal subtraction renormalization scheme, renormalization 
is accomplished by subtracting the pole terms from the coefficients and then taking the limit 
e — > 0. The functions of u in (|62|) and ( |63|) can be expressed in the form 



f(u,e) = J e ■ , 

u + e 



(64) 



where F(u + e, e) is analytic in u + e at u + e = 0. The renormalization of this expression 
by minimal subtraction involves subtracting the poles in e from the coefficients in its power 
series expansion in u, followed by taking the limit e — > 0. The final result is rather simple: 



/ats(u) 



u 



F(u, 0)-G(u) 



where G(u) is the Borel transform of the function 

(47r) x r(4 + 2x) 



G{x) 



6r(l -x)r 2 (2 + a; 



-F(0,-x) . 



(65) 



(66) 



We can derive this result through several simple steps that allow us to take the limit e —>■ 
in parts of the expression (|B~4]) . The first step is to subtract e^ A ~ c ^ e F(0, e)/(u' + e) from 
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e ( A c ^ u ' F(u' + e, e)/(u' + e) in fl64l) , and add it back in again. The subtraction cancels the 
pole at u' + e = 0, so we can take the limit e — > except in the added term: 



f(u,e) 



F(u,0) - e u ( A - c )F(0, 0) 



■it 



+ F(0, e) e - u/e+ ^- c)t exp | 



tt 



ff( e ) „e(d/8u') 



d_ 

du' 



it' 



It is convenient to use the following integral representation for l/(u' + e) : 



1 



u> 



dae- {u ' +e)a 



(67) 



(68) 



Inside the integral, the action of the differential operator d/du' is simply multiplication by 
—a. After setting u' = u, the action of the exponentiated operator in flBTD can be carried 
out analytically: 

d 



exp 



u m e <9/au') _ " 
e du' 



1 



u' + e 



1 

ug(e) 



(69) 



Inserting ( |5DD into flB7|), we can take the limit e — > in the term with the factor e ug ^l e . The 
resulting expression for the Borel transform is 



f(u,e) 



1 

tt 



F(«,0) - ^lil e -A + (A-c), 



Minimal subtraction can now be carried out by using the identity 



G(-e)e 



-u/e 



MS 



G(u), 



(70) 



(71) 



where G(u) is the Borel transform of the function G(x). Inserting this into (70), we obtain 
our final result (|65]). 

We now apply (^) to the expression (^) to obtain an expression for the Borel transform 
of the pole mass in terms of the running mass m(/i) defined by minimal subtraction: 

' — 1 II i 1 ; / r I I . 1 ■ .' / r 

G m {u) 



m 



poie(*) = m(ji) <6(t) + 



nu 



r(i + «)r(3 - 2u) 



e c m 2 (/j,) J (l-2u)T(3 



u 



where G m (u) is the Borel transform of the function 

(3 + 2x)r(4 + 2x) 



G m (x) 



(72) 



(73) 



9r(l-a;)r 2 (2 + x)r(3 + x) ' 

This agrees with Eqs. (4.3) and (4.4) of Ref. 0. The expression (|72|) for m po i e has poles 
in the Borel parameter t = u/(3q. These poles are called renormalons. They are classified 
as ultraviolet renormalons or infrared renormalons according to whether they arise from the 
integration region of high loop momentum or low loop momentum. The pole in (|72|) from 
high loop momentum are those that appear in the factor T(l+u). Thus there are ultraviolet 
renormalons in the pole mass at the negative integers u — — 1, —2, • ■ -. There are infrared 
renormalons at the positive half-integers u = |, f, - • •, and also at u — 2. From (|BB"D, we 
see that SZ has renormalons at the same positions as the pole mass except that there is no 
renormalon at u — — 1. 
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C. Quarkonium-to- Vacuum Matrix Elements 



Now we calculate the Borel transforms of the short-distance coefficients in the NRQCD 
factorization formulas for the quarkonium-to-vacuum matrix elements. We define the short- 
distance coefficients C r (a s ), where T = V, A, T, P, by the factorization formulas 



(0|*7^|J/^) 

(0\^a^\J/ijj) 
(Ol^^c) 



s* [i + c v (a s )} (o| x Vvi JM + ••• , 

S*> [l + C A (a s )} (0\ X ^\Vc) + ■•• , 
-i{S"° 5™ - 8 v0 8*) [l + C T {a s )] (0|xVVlW> + 
1 + C p (a s )} (0|xVl^c) + ••• • 



(74) 
(75) 
(76) 
(77) 



We use T to represent the matrices 7 M , 7^75, cr^, and 75, so that the matrix elements on 
the left side of (f74D-(f77|) are denoted genericly by (0|\[>r^|F). 

The short-distance coefficients C r (a s ) in (f74|)-(|7"7|) can be determined by matching the 
corresponding equations for |cc) states computed in perturbation theory. The simplest 
choice is a state |cc) that consist of a c and c at rest. We must compute (Ol^r^lcc) using 
perturbative QCD and (0|%^ cr l ip\cc) and (0|xV|cc) using perturbative NRQCD, and then 
determine C r (a s ) by matching the left and right sides of (|74|)-(|77|). By Borel transforming 
both sides of the equation before carrying out the matching, we can determine the Borel 
transform C r (t). 

We first consider the QCD side of the matching conditions. The c and c at rest have 
equal four-momenta p = (m, 0) and are represented by the spinors 



u(p) 




v(p) 




(78) 



The diagrams that dominate the radiative corrections in the large- Nf limit are obtained by 
inserting a chain of light quark bubbles into the gluon propagator in the one-loop diagrams. 
The diagrams with propagator corrections on the c and c lines reduce to the tree level 
matrix element vTu multiplied by 5Z(a s ). After Borel transforming, they give 5Z(t)vTu. 
The remaining diagrams are vertex corrections, such as the diagram for cc — > e + e~ in Figure 
4. The Borel transform of the vertex correction is 



-P 



,1+u 



v(p) 7 1 



1 



1 



/— i> — m + it /+ ^ — m + ie 



7„«(p). (79) 



The dimensionally regularized integral in (|79|) can be evaluated by introducing Feynman pa- 
rameters. After a straightforward calculation, we find that (|79| ) reduces to D r (t) v(p)Tu(p), 
where 



3 ( 7 -5/3)e 



37T 



Jc 



e c m 2 



" T(u + e)T(l -2u-2e) 

(1 + 2u + 2e) r(3 - u - 2e) 
e)(3 + 2e) + {l-2e)u} 



(l_2c-«) [(1 
(3-2e) [(1 -e) + m + u 2 } 
2(1 - e)(l - 2e) - (3 - 3e + 2e 2 )u + 
(3-2e)(l + e + u)(2-2e-u) 



2e)u 2 



r 
r 
r 
r 



v 

A 
T 
P 



(80) 
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where the operator J e acts on every function of u to the right. The prescription we have 
used for 75 in D dimensions is that {75, 7 M } = for all \i. 

We next consider the NRQCD sides of the matching conditions. If we use dimensional 
regularization, the only scale in the perturbative NRQCD matrix elements is the relative 
momentum q. Having set q = 0, the only matrix elements that survive are (O|xV"0|cc) and 
(0|x^|cc). The matrix elements of all higher dimension operators vanish. The radiative 
corrections to the matrix elements vanish in dimensional regularization, since there is no 
momentum scale in the loop integral. Thus the NRQCD matrix elements are trivial: 



(0|x f <xV|cc) = 2m7?V£ 
(0|x^|cc) = 2mry^ . 



(81) 
(82) 



Since these matrix elements are independent of a s , the Borel transforms of the right sides 
of (|T4])-(|T7|) are obtained simply by replacing 1 + C r (a s ) by 5(t) + C v (t). 

By matching the Borel transforms of the QCD and NRQCD sides of ( 
that the short-distance coefficients are given by 

C r (t) = D T (t) + 5Z{t) , 



-(|77|), we find 



(83) 

where 5Z{t) is given in (|63j) and D T (t) for T = V, A, T, and P is given in (|80|). The 
Borel transforms C v and C A are finite at u + e = and therefore do not need to be 
renormalized. This is simply a consequence of the conservation of the vector current and 
the partial conservation of the axial current. We can therefore set e = in these coefficients: 



C v (t) 
C A (t) 



2 


( " 2 


3tt 


\e c m 






7T \ 


e C m 2 



4 + 3u r(l + u) T(3 - 2u) 



(l-2u)(l + 2u) r(3-tt) 

1-u-u 2 r(i + «)r(i - 2u) 



l + 2u 



r(3 



u 



(84) 
(85) 



In the Borel transforms C T and C p , there is a pole at u + e = 0. This pole indicates that 
there are ultraviolet divergences in the coefficients of the power series expansion in u. These 
divergences are removed by multiplicative renormalizations of the QCD operators ^a^^ 
and ^75^/. If we use the minimal subtraction renormalization scheme, the renormalization 
can be carried out by making the substitution ( |65D [[ CHANGE EQ. NO. d65|). ]] Our final 
results for these coefficients are 



C T (t) 
C p (t) 



1 



3iru 



i + iom + 6m 2 r(i + «)r(3-2«) 



e c m 2 



2 

7TU 



I (l-2u)(l + 2u) T(3-n) 
l-2u + 2u 2 T(2 + u) T(l - 2u) 



G T {u) 



e c m 2 



1 + 2u 



r(3 



u 



(86) 
(87) 



where G T {u) is the Borel transform of the function 



G T (x) 



(1 + 2x)r(4 + 2x) 



3r(i-x) r 2 (2 + x)r(3 + x) 

and G m (u) is the Borel transform of the function G m (x) defined in (|73| ). 



(88) 
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From (|84D -(p7D, we see that all the coefficients have ultraviolet renormalons at the neg- 
ative integers u = — 1, —2, • • •, except for C p , which does not have a renormalon at u = — 1. 
They all have infrared renormalons at the positive half-integers u = |, |, • • •, and at u = 2. 
The coefficients C A and C p have an additional infrared renormalon at u = 1. Finally, all 
the coefficients have a renormalon at u = — |. This appears to be an infrared renormalon 
in spite of the fact that it appears on the negative real axis, because it arises from the 
low momentum region of the integral over the loop momentum I in (ff9|). However the low 
momentum region of that integral is actually canceled by the low momentum region of an 
integral on the NRQCD side of the matching condition. The integral on the NRQCD side 
has the form 

16nm ( /i 2 V r dP~ 2 H 1 

3" [e*) J (2vr) 3 - 2e Jpf+^ ' (89) 

This vanishes in dimensional regularization due to a cancellation between an infrared pole at 
u + e = — \ and an ultraviolet pole at u + e = — |. Since the infrared pole is canceled by the 
QCD contribution to C r , the renormalon at u — —\ is actually an ultraviolet renormalon 
of NRQCD. 

To obtain the power series expansions for the coefficients C r (a s ), we expand the coef- 
ficients C r (t) as a power series in t and then invert the Borel transform by substituting 
t n ~ l — > (n — The first two terms in the expansion of C v (a s ) are 

Wi \ _ 8 oisifi) , /_ 16 fi \ (3 Q a 2 s 



C v (a s ) = 2__li og ^ p^i. (90) 

The coefficient of a s was first calculated by Barbieri et al. | 20|| . From the a 2 term in (p0|), 



we can read off the Brodsky-Lepage-Mackenzie (BLM) scale [^l[j for this process, which is 
the scale /x for which the correction proportional to Nfa 2 vanishes. We find /i = 1.46 e c / 2 m. 



D. Decay of r\ c into 77 

We next calculate the Borel transform of the short-distance coefficient of (0|x^|^ c ) in 
the NRQCD factorization formula for the T-matrix element for 7] c — > 77. The factorization 
formula with tree-level coefficients is given in ([32]) . Up to this point, it has not been necessary 
to define carefully the mass m that appears in the factorization formula. It could equally 
well be the running mass m(/i) defined by a minimal subtraction renormalization scheme, 
or the pole mass m po i e defined by the location of the pole in the propagator in perturbation 
theory, or any other mass that differs from these by perturbative corrections. However, in 
order to calculate the short- distance coefficient beyond tree level, it is necessary to specify 
the definition of m. The choice is somewhat arbitrary, since a different choice can be 
compensated by changes in the perturbative corrections to the short-distance coefficient. 
We choose to use the pole mass. We define the short- distance coefficient C 77 (a s ) by the 
factorization formula 

V-n = ~i ,TlT n 2 x e 2 • r [ 1 + C^(a s ) } (0\ X ^\Vc) + ■ ■ ■ , (91) 
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Our choice of the pole mass is dictated primarily by the algebraic simplicity of the final 
result for the Borel transform of C 77 (q; s ). The coefficient C 77 in (plf) can be determined 
by matching the perturbative T-matrix elements for cc — > 77 in QCD and in NRQCD, as 
in (H). If we take the c and c to be at rest, their invariant mass is M cS = 2m pole , and the 
matching condition reduces to 



- 1 cc— >77 



= ei x e 2 • f [ 1 + C^(a s ) } (0\ X ^\cc) 

pQCD ^polc 



+•••. (92) 

p NRQCD 



By Borel transforming both sides of the equation, we can determine C 77 (£). 

We first consider the NRQCD side of the matching condition. Having chosen the relative 
momentum of the cc pair to be 0, the radiative corrections to the matrix element vanish in 
dimensional regularization. The matrix element therefore reduces to 

(0|xV|cc) = 2m pole r^ . (93) 

The Borel transform of the right side of (|9"2"|) is therefore obtained by replacing 1 + C 77 (a! s ) 



by 5(t)+C^(t). 

We next consider the QCD side of the matching equation fl9"2"|). Since the c and c are 



at rest, they have equal 4- momentum p = (m, 0) and their spinors are given in ([781). The 
diagrams for the T-matrix element in the large Nf limit are those given in Figs. 5a-e, 
together with the diagrams with self-energy corrections on the external c and c lines. In 
each of the diagrams, the spinor factor can be reduced to r\vYt v u^ where r = (k% — &2)/2. 
The contribution to C 77 (t) can be obtained from the Borel transform of the diagram by 
removing the factor 

Q 2 e 2 

——^ex^urxv^Y^^p) = -iQ 2 e 2 e 1 x e 2 • . (94) 

The ultraviolet divergences associated with vertex corrections and wavefunction renor- 
malizations cancel by the Ward identities of QED. The only renormalization that is necessary 
is therefore mass renormalization. It is convenient to use on-shell renormalization. Explicit 
regularization of the loop integrals is unnecessary, since the Borel parameter itself provides 
a regulator. The mass counterterm and the wavefunction renormalization constant can be 
obtained from fl62|) and ( |63"|) by setting e = 0: 



^ = m £(^-lz!iSlzWt20, (95 ) 

7i \e u m z J u i (3 — u) 

,z w = - 1 {4-X r(i - r y r V" ) . (96) 

7i ye°77r J u 1(6 — u) 

The short-distance coefficient C 77 (t) can be written 

C^it) = C^t) + C 2 (t) + C 3 (t) + 5Z(t) , (97) 

where C\(t) is the contribution from the box diagrams in Fig. 5a, (^(t) is the contribution 
from the vertex corrections in Figs. 5b and 5c, Cs(t) is the sum of the contribution from 
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the propagator correction in Fig. 5d and the mass counterterm in Fig 5e, and 8Z(t) is the 
contribution from propagator corrections on the external lines. 

The sum of the Borel transform of the box diagram in Fig. 5a and the diagram obtained 
by interchanging the photon lines is given by 



2„2 



16nQ 2 e 



x 

The integral in Eq. 
(|94|), we obtain 

Ci(t) = 



e c J Ji (-py+ u 



7 +7 



1 



/+ / — m 



/— / — m 



/— f> — m 



r 



/+ p" — m 



lau{p) 



(98) 



can be carried out using Feynman parameters. Removing the factor 
" (2 - u) (1 + 2u + 3m 2 ) T(l - 2u)r(l + u) 

(99) 



3n \ e c m 2 



u(1 + u)(1 + 2m) 

i r(i- u )r(i±M) 



r(3 - u) 



i+u r ( ^ 



The Borel transform of the vertex correction in Fig. 5b is obtained by replacing 7 M by 



. 167T / fl 



2\ u 



:,C 



7 



r 



i (-Z 2 )^" ' /- /-m ' /+ 



7a 



(100) 



Using Feynman parameters to evaluate the integrals for the four vertex correction diagrams, 
we find that their contribution to C 77 is 



C 2 (t) = - 



37r V e c m 2 



2(1 -2m) r(i-2it)r(i + «) 



ii- 



r(3-«) 



+ - 

u 



i r(i - «)r i r(2-«)r(a±a 



r (¥) 



2T 



4— M 

2 



(101) 



The Borel transform of the self-energy correction for the internal quark line in Fig. 5d is 
given by (|60D with p replaced by — r. With a straightforward calculation, we obtain 



1 



E/(-/,t) = 



3/T V e c 'm 2 



rf-u^r 



2+M 

2 



r(i=H) 



+ 2m 



1+M 

2 



r(¥) 



(102) 



The diagram in Fig. 5e has the mass counterterm Sfh inserted into the internal quark line. 
Adding the contributions to C 77 from Figs. 5d and 5e, we obtain 



Cs(t) = - 



37r \ e c m 2 



1 2r(i - w )r 3(1 -u) r (i - 2u) r(i 



F / 3— » 
1 1 2 



U 



r (3 - u ) 



. (103) 



Adding up (p9f) , (|101|) , ( |103| ), and (|9~6D, the Borel transform of the short-distance coefficient 
for the process r] c — >• 77 is 
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Crr(t) = — 

D7T 



e c m 2 



2(1 + 2u - Au 2 - 5u 3 + 3u 5 ) T(l - 2u)r(l + u) 



u 2 (l + u)(l + 2u) r(3-u) 

i r(i- M )r(^) i r(i-n)r(^) 



1 + it 



r(s=*) 



r(¥) 



(104) 



While individual terms in ( |104| ) have a double pole at u — 0, there are cancellations that 
make the expression analytic at u = 0. From eq. ( |104|) , we see that there are ultraviolet 
renormalons at the negative integers u = — 1, —2, • • •. There are infrared renormalons at 
the even and odd positive half-integers u = |, 1, |, 2, • • •. Finally, there is the NRQCD 
ultraviolet renormalon at it = — ~. 

we expand ( |104j ) in powers of u and 



To obtain the power series expansion of C 



TO ( 



carry out the inverse Borel transform by substituting t n 1 
in the expansion are 
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n — 1)! a™. The first two terms 
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(105) 



The coefficient of a s was first calculated by Barbieri et al. [p2]| . From the vanishing of the 
coefficient of /3o« 2 , we identify the BLM scale for this process to be fx = 4.49 e c ' 2 m. 



IV. RENORMALON SINGULARITIES IN QUARKONIUM ANNIHILATION 

DECAYS 

We have calculated the Borel transforms of the short- distance coefficients C r (a s ) for 
T = V,A,T,P, and 77 in the large- Nf limit. They have renormalon singularities in the 
form of poles on the real axis of the Borel parameter t = u//3q. There are ultraviolet 
renormalons at the points u = —1, —2, • • -, infrared renormalons at |, |, 2, |, 3, • • *, 
and a NRQCD ultraviolet renormalon at u = — |. Each renormalon generates n! growth 
of the coefficients in the perturbation series, as shown in (|48|). The asymptotic behavior 
of the coefficients is dominated by the renormalons closest to the origin at u — ±|. The 
contributions from a higher renormalon at u = k/2 are suppressed by k~ n . In this section, 
we study in detail the dominant renormalons at u = +\ and u = — |. We show that 
the Borel-resummable renormalon at u = — | arises from the cancellation of the Coulomb 
singularity in the short-distance coefficients and is universal in the sense that its residue is 
identical for all color-singlet S'-wave decay processes. We show that the ambiguities from the 
renormalon at u — % can be absorbed into the nonperturbative NRQCD matrix elements 
that contribute through relative order v 2 . We also deduce the residues of the renormalons 
at u = ±| for the decay J/ip — > 777. 

A. Universality of the u = — ^ renormalon 

The Borel transforms of the short-distance coefficients C r (a s ) are given in (|g^])-(|gT|) for 
r = V, A, T, P and in (|104| ) for T = 77. By examining these expressions, we see that they 
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all have a renormalon at u = — 5 and that the residues of the poles are identical for all T. 
Thus they can all be written in the form 

^rrW(-^?H rw > (106) 

where B r (t) is analytic at /3ot — — |. Since the pole at u = — \ is located on the negative 
real axis, the contribution of this renormalon to C r (a s ) is Borel resummable. If B r (a s ) is 
the inverse Borel transform of B r (t), then the inverse Borel transform of (|106|) is 

C T (a s ) = -f a s / dx — f- + B\a s ) . (107) 

in /i jo 1 + 2p a s x 

We have suppressed the /i-dependence of the coupling constant a s (n). For sufficiently small 
values of a s , the integral is approximately 1. However its power series expansion in a s is an 
alternating series whose coefficients diverge like n\. 

We found that the residues of the renormalon at u — — \ are the same for all processes 
that we calculated. We will show that this result is general, that the residue is identical for 
all color-singlet S"-wave annihilation decay processes. The appearance of this renormalon 
can be traced back to a remnant of the cancellation of the Coulomb singularity in the 
short- distance coefficients. 

Any S'-wave charmonium decay proceeds through the annihilation of the cc pair. The cc 
annihilation amplitude at tree level defines a vertex T(p,p), which depends on the 3-momenta 
of the c and c. The scale of the momentum dependence is m. The vertex can therefore be 
approximated by a momentum independent vertex when c and c have 3-momenta that are 
small compared to m. The QCD radiative corrections to the annihilation amplitude include 
the diagram in Fig. 4 in which a virtual gluon is exchanged between the c and c. If the 
relative momentum between the c and c is zero, there is a pinch singularity in the loop 
integral. This is the Coulomb singularity. If the c and c have a small relative momentum 
mv, the singularity appears as a 1/v term in the annihilation amplitudes. Since NRQCD by 
construction matches full QCD at momenta small compared to m, this Coulomb singularity is 
the same in the two theories. It therefore cancels in the short-distance coefficients. However, 
the cancellation leaves a finite remainder. If we insert a chain of n bubble diagrams, into the 
propagator of the exchanged gluons, there is still a cancellation of the Coulomb singularity, 
but the remainder grows like n\. These remainders form an alternating series, which is Borel 
resummable. This series produces the renormalon in the Borel transform at u — — \. 

We proceed to calculate the Borel transform of the sum of diagrams in Fig. 4 in the 
full theory. The expression is identical to (^), except that T is replaced by a momentum- 
dependent vertex T(p — l,p + I): 




+ m T-T(p-l,p + l)T° J + J +m . 7 „ 



u{p) . (108) 



P-21-p + ie ^ ,r ' l 2 + 2l-p + ie 

Since the Coulomb singularity comes from the region of I <C m, the / terms in the numerator 
and in the vertex T can be neglected. Using phi = mu and v p 1 = —mv, the above equation 
can be reduced to 
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(2 \ u 
v{p)T a T(p,p)T a u{p) 

dog) 

i (-1 2 - ie) 1+u I 2 - 2ml + ie I 2 + 2ml + ie ' y ' 

If the annihilation vertex T is a color- singlet, the factor vT a TT a u reduces to (4/3)vTu. If V 
is a color-octet vertex, the factor 4/3 is replaced by —1/6. 

To evaluate this integral, we first integrate over l Q by using the contour integral method. 
Closing the contour in the upper-half plane, the Coulomb singularity appears in the contri- 
bution from the pole at 1$ = —m + y/ m 2 + l 2 + ie. Keeping only this term, (|109|) reduces in 
the case of a color-singlet annihilation vertex to 

^ (£) m rfep) u(p) / wf { 2 m v^-2^r ' (110) 

The first factor in the integrand reduces to l/(l 2 ) 2+u for small 1. It is easy to see that if 
( P-lOp is expanded as a power series in u, all the coefficients have a linear infrared divergence 
coming from the small 1 region. This just corresponds to the Coulomb singularity. 

We next consider the Borel transform of the sum of the diagrams in Fig. 4 in NRQCD. 
Since there are no momentum scales in the diagrams and since it must match ( |110| ) at low 
momenta, the result must be 

(^) v(p) r(p, P ) u{p) J ^-L ^ . (in) 

The contribution of the Coulomb singularity to the short- distance coefficient C r (t) is ob- 
tained by subtracting (|111|) from ( |110| ) and removing the factor vTu. The Coulomb singu- 



larity cancels, so the power series expansion of C v {t) has finite coefficients. However, the 
integral in QUI]) has a power ultraviolet divergence when u > — ^, and this produces a pole 
in C T {t) at u = — |. The pole can be isolated by imposing a cutoff |1| > tz on the integral: 

ff(t) . J*= g)" L |L . (112) 

Evaluating the integral, we find that the residue of the pole is identical to that in ( |106| ). 

It is easy to see that the renormalon at u = — \ and its residue are gauge invariant. If we 
had used the gluon propagator (|55|) to carry out the calculation in a general covariant gauge, 
there would be additional terms in ( |108|) with 7^ replaced by J/. Making the decomposition 
f — (p 1 — m) — (p 1 — )/ — m), the first term vanishes due to the equation of motion for u(p), 
while the second term cancels the quark propagator. This removes the pinch singularity, 
and therefore these terms have no pole at u — — 5- 

The residue of the u = — \ renormalon is universal for all color-singlet S'-wave quarkonium 
annihilation decays. The universality follows from the fact that in the low-momentum region 
of the loop integral for the QCD amplitude (|108|), the cc annihilation amplitude can be 



approximated by a vertex T(p,p) that is independent of the loop momentum and nonzero 
for S'-wave annihilation. The low-momentum regions of the QCD loop integral and the 
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corresponding NRQCD loop integral must match. But there are no momentum scales in 
the NRQCD loop integral, and therefore its ultraviolet behavior is also fixed. The residue 
of the ultraviolet NRQCD renormalon at u = — | depends only on whether the annihilation 
amplitude T(p,p) is color-singlet or color-octet. Note that this renormalon is generated only 
by loop diagrams in which a gluon is exchanged between the c and c. Other loop diagrams, 
such as quark self-energy diagrams and the diagrams in which a gluon connects the core 
line to the annihilation vertex, give no contributions to the renormalon simply because there 
is no pinch singularity in the QCD diagrams. 

B. Cancellation of the u = i renormalon ambiguities 

The Borel transforms C r (t) in (|34"D-(|S7|) and ( |104j ) have an infrared renormalon at u = 
+|. They can therefore be written in the form 

^ ^ r^kt + ir « • (113) 

where A r (t) is analytic at Pot = | and the residues of the poles are 



9tt e c ' 2 



m 



a T = — , (116) 



3tt e°l 2 
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or = --<£-, (H7) 

a 9ire c /*m- U * } 

Since the pole at u = | is located on the positive real axis, the inverse Borel transform is 
not unique. Taking the principal value of the integral in ( |47| ) is only one of the possible 
prescriptions. If A r (a s ) is the inverse Borel transform of A r (t), the general expression for 
the inverse Borel transform of ( |113| ) is 

C r (a s ) = a r a s p/°° dx + A r (a s ) - K^fe- 1 /^ , (119) 

Jo 1 — 2p a s x p 

where P/ denotes the principal value integral, K is an arbitrary constant, and a s = a s (n). 
The last term in ( |119| ) represents the ambiguity in the inverse Borel transform, which is 
proportional to the residue of the pole in the integrand. Using the one-loop expression (|50|) 
for the running coupling constant a s (fi) and the fact that aF includes a factor of yu/m, we 
see that the ambiguity in ( |119| ) is proportional to Agcc/ m - This ambiguity arises because 



summing the series of loop diagrams with bubble chains inserted into the gluon propagator 
is equivalent to replacing the strong coupling constant by a running coupling constant at 
the scale of the gluon momentum. However, the perturbative running coupling constant has 
a Landau pole at the momentum scale Aqcd- This pole appears in the integration region 
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and therefore gives rise to an ambiguity in the loop integral. The renormalon ambiguity 
implies that the short-distance coefficients defined by minimal subtraction necessarily have 
some sensitivity to nonperturbative effects involving the energy scale Aqcd- 

The infrared renormalons give ambiguities in the short- distance coefficients in NRQCD 
factorization formulas. The decay rates for J/ip — > e + e~ and r] c — > 77 and the matrix 



elements in (|Tq)- (|77D are physical quantities, so they cannot have any ambiguities. The 
NRQCD factorization formulas simply represent the separation of scales. All effects of mo- 
mentum scales m and larger are included in the short-distance coefficients, while effects of 
lower scales, such as mv, mv 2 , and Aqcd, are absorbed into the NRQCD matrix elements. 
Since infrared renormalons arise from small momentum scales, any ambiguities in the coef- 
ficients due to infrared renormalons must be compensated by corresponding ambiguities in 
the nonperturbative matrix elements. Thus it must be possible to consistently absorb the 
ambiguities proportional to K in (|119|) for V = V , A, T , P and 77 into various NRQCD 
matrix elements. 

The ambiguities in a coefficient from an infrared renormalon at u = k/2, where k is an 
integer, is proportional to (Aqcd/W) & . Since the NRQCD matrix elements involve several 
different scales, it is not immediately obvious which matrix element an ambiguity propor- 
tional to Aq CD should be absorbed into. We will make a simple hypothesis that is consistent 
with the velocity-scaling rules of NRQCD and verify that the ambiguities due to the u — +| 
renormalon are consistent with this hypothesis. Our hypothesis is that ambiguities propor- 
tional to Aq CD can be absorbed into NRQCD matrix elements that are suppressed by v 2k 
or less. If our hypothesis is correct, ambiguities due to the u = +| renormalon should be 
absorbed into the leading order matrix elements and those suppressed by v 2 . In the case of 
the J/i/j, these matrix elements are (0\x^cr l ip\J/?P) and (0|x^D 2 cr l '?/'| J /if)) ■ In the case of rj c , 
they are (0|xVl^c) and {0\x^~D 2 ip\ri c ) . Because of the spin-symmetry relation (|7|), there are 
only three independent matrix elements that are suppressed by v 2 or less. It is convenient 
to introduce the following shorthand notation for these three independent matrix elements: 

(jq,) = e -(OlxWl JMe)>, (120) 
(JC V ) EE (0\ X ^\Vc), (121) 

(£i»> = — 2 (0|x f DVl»7e> • (122) 
m z 

We first consider the cancellation of the renormalon ambiguities in the QCD matrix 
elements in ([73])- (177]). The coefficients of the matrix element {JCjj 2 ) i n these factorization 
formulas are given in ( |14|) and ((42j)- fl44]) . The coefficients C T are functions of /3oa s divided 
by /3q and are therefore suppressed by 1/Nf in the large- Nf limit. Since the ambiguities AC r 
in the coefficients are suppressed by l/Nf, the ambiguities A(fC^), A(JC V ), and A(JCo2) in 
the matrix elements must also be suppressed by l/Nf. The cancellation conditions for the 
renormalon ambiguities at leading order in l/Nf are therefore linear relations: 

AC v (JCtp) + A(/Q> + ^A()C D2 ) = 0, (123) 
AC A ()C V ) + A(IC V ) + ^A(/C D2 ) = 0, (124) 
AC T {1C^) + A(/C V) ) + 1 -A(1C D 2) = 0, (125) 
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AC P (JC V ) + A{)C V ) = . 



(126) 



From ( 119 ), the ambiguities in the coefficients are 

7ia T A 



AC 



-K 



QCD 



. (127) 

Po A 4 

with the coefficients aF given by (|114|) -( 117 ). The cancellation conditions ( 123|) - fll26|) are 
four equations for the three unknown ambiguities in the matrix elements, which provides a 
consistency check. The solutions are 



A(/C„> 
A(/C D2 ) 



—K 
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QCD 
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QCD 



(128) 
(129) 
(130) 



3 /3 e c l 2 m 

These results can be verified by perturbative calculations of the large- Nf limits of the Borel 
transforms of the NRQCD matrix elements. 

Now we examine the renormalon cancellations for the process r\ c — > 77. From fl32"D and 
(PT|) , the factorization formula for the T-matrix element can be written in the form 



2Q 2 e 2 
-i — — ei x 62 ■ r 



M, 



Vc 



1 + C 



77 



(131) 



We have used the identity (0) to eliminate the pole mass m pole in the prefactor in favor of 
the meson mass M Vo . The cancellation condition for the ambiguities from the renormalon 



at u 



is 



AC^{JC V ) + A(1C V ) + -A(K D 2) = 0. 



(132) 



Using the value of a 77 in ( |118| ) and the values of A(JC V ) and A{1C D 2) given in ( |129| ) and 
( P-30Q , we see that ( |132j ) is satisfied. 

We can use the results given above to verify the cancellation of renormalon ambiguities in 
the Gremm-Kapustin identity, which is given at tree level in (|T7|). The correct generalization 
to higher order in a s is to replace m by m po i e . The analogous identity for the matrix elements 
of rjc is 



(0\x^\Vc) 



-m 



pole 



(M, 



Vc 



2m 



pole. 



(0|xVhc) (l + 0(v 



(133) 



From the expression ( [72] ) for the Borel transform of m po i e , we find that the ambiguity due 
to the renormalon at u 



+5 is 
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pole 
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QCD 



(134) 



The leading contribution to the ambiguity on the right side of ( |133j ) comes from the factor 
M nc — 2m po i e . The ambiguities from the factors on the right side of (|133|) are suppressed by 



v 2 from the factor M Vc — 2m po \ e . Thus the renormalon cancellation condition is 



A(/C D2 ) 



2 Ampoie^ 



m 



(135) 



Using (|130Q and ( |134| ), we see that this is satisfied. 
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C. Leading Renormalons for J/ip — ► 777 
The one-loop QCD correction to the short-distance coefficient for the process J/tp — > 777 



has been calculated by Mackenzie and Lepage p3[ . The correction to the decay rate is a 
multiplicative factor 1 — I2.6a s /ir, which gives a correction of about —140% if we take 
the coupling constant to be a s (m) ~ 0.35. Since it is the leading order correction, the 
coefficient of a s is independent of the renormalization scale and the renormalization scheme. 
The enormous radiative correction raises questions about the origin of the large coefficient 
of a s and about the convergence of the perturbation series. Some insights can be gained 
by studying the renormalon singularities for this process. As shown above, the large order 
behavior of the short- distance coefficient is governed by the u — ±| renormalons. We know 
the residue of the u — — | renormalon in the large Nf limit, since it is universal for all S- 
wave annihilation decay processes. The residue of the u — ^ renormalon can be determined 
by requiring the cancellation of the renormalon ambiguities in the decay width. Knowing 
the residues of the leading renormalons, we can estimate the asymptotic behavior of the 
perturbation series. 

The decay width for J/ip — > 777, with the coefficient of the first relativistic correction 
calculated to leading order, is 



r(j/i/> -> 777) 



16(tt 2 - 9)Q 6 a 3 



1 + C^(a s )) (0|*W|J/V> 



(136) 



24(vr 2 - 9) m? 

The coefficient of the matrix element (0|x^D 2 en/>| J/ip) was first calculated by Keung and 



Muzinich in numerical form [24] and by Schuler in analytic form [25]. The condition for the 
cancellation of the u = \ renormalon ambiguities in Y^J/ip — > 777) reads: 

771-2 _ 94 

AC 3 ^> + A(X^> + 24(?r2 _ 9) A(/C D2 ) = 0. (137) 

Using the values of A(/C^) and A(/C/j2) determined in ( \L28j) and (|130|) , we obtain 

AC^ = K U f ~ 8 l 4^ ■ (138) 
9(vr 2 - 9) (3 e c / 2 m K ' 

Comparing with (|127| ), we can read off the residue a 37 of the u = \ renormalon, which is 
defined by ([1131): 

This residue is rather large. It differs from the residues a v and a 77 for the decays J/if> — > e + e~ 
and t] c — > 77 by factors of 5.4 and —11.8, respectively. The large one-loop radiative correction 
may be related to the large residue of the u = \ renormalon. 
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V. ESTIMATES OF HIGHER ORDER RADIATIVE CORRECTIONS 



One of the problems encountered in phenomenological applications of the NRQCD fac- 
torization formalism is that the one-loop corrections to the short- distance coefficients for 
several decay rates are uncomfortably large. This raises questions about the convergence 
of the QCD perturbation series. The calculations of the short-distance coefficients in the 
large- Nf limit can shed some light on the origin of the large radiative corrections. They can 
also be used to provide simple estimates of higher order radiative corrections. It has been 
found empirically that if we take the coefficients of the perturbation series in the large- Nf 
limit and make the shift Nf — > Nf — 33/2, we obtain reasonable approximations to the higher 
order coefficients in cases where they are known exactly This prescription, which has 
been called "naive nonabelianization" , puts the Borel singularities at the correct points on 
the t-axis. We will use this prescription to give simple estimates of higher order corrections 
in the perturbation series for C v (a s ), C 77 (a s ), and C 3l (a s ). 

The Borel transform C r (t) for the short- distance coefficients in the decay rate for J ftp — > 
e + e~ and rj c — > 77 are given in ( |84] ) and ( 104| ). To obtain the perturbation expansion for 
C r (a s ) in the large- Nf limit, we expand C r (t) as a power series in u, and then invert the 
Borel transform by substituting t n ~ l — ► (n — 1)! a". We set j3 = (33 — 2Nf)/12n and take 
Nf = 3 for charmonium decays. We choose the MS scheme by setting C = —5/3 and we 
take the scale of the running coupling constant to be \i = m. The resulting perturbation 
series are 



C v {a s ) 



-2.7a - 5.5 a 2 - 150 a 3 - 870 a 4 - 38000 a 5 - 360000 a 6 + 
-1.7 a + 5.1a 2 - 25 a 3 + 1100 a 4 - 900 a 5 + 520000 a 6 + ■■ 



(140) 
(141) 



where a = a s (m)/n. We have kept only two significant figures in each of the coefficients. 
The behavior of the perturbation series (|140|) and (|141| ) can be understood by examining 



the contributions from the renormalon singularities at u = | and u = —-k. Asymptotically, 
these renormalons dominate the coefficients in the perturbation series, because the contribu- 
tion to the coefficient of a™ from a higher renormalon at u = k/2 is suppressed by k~ n . The 
contributions to the perturbation series ( |140p and (141) from the renormalons at u 
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(n-l)![2/3 a s (/i)] n , (142) 
n-l)\[2f3 a s (fi)} n . (143) 



I renormalon is 



( |143| ) are proportional to the residues of the u 
of the two renormalons depends on the renormalization scale \i. The u 



The alternating series in (|142|) and (|143|) are identical, since the u = 
universal for all S"-wave annihilation decay processes. The same-sign series in (|142|) and 

+| renormalon. The relative importance 

— \ renormalon 

dominates if fi m and the u — +| renormalon dominates if w ^> m. The choice of fi 
that gives equal contributions from the two renormalons is [i = lA8e c ^ 2 m for (|142|) and 
H = 2.19e c ' /2 m for ( |143| ). In the MS scheme, these scale are 0.64m and 0.95m. These values 
suggest that fi = m is a reasonable value for the renormalization scale in the MS scheme. 
For comparison, the corresponding BLM scales are /1 = 0.63m and \x = 1.95m, respectively. 



28 



In the MS scheme (C = — |) with the scale \i = m, the first few terms in the series (|142|) 
and ( |143| ) are 



C v (a s ) 4.0a - 7.4a 2 - 160a 3 - 900a 4 - 39000a 5 - 370000a 6 + ••• , (144) 

C 77 (a s ) ~ 0.12a + 11a 2 + 4.8a 3 + 1300a 4 + 1200a 5 + 540000a 6 + ■•■ . (145) 



Comparing (|140|) and ( 144 ), we see that the leading renormalons begin to dominate the 



large- Nf coefficients of C v at the a 3 term. Comparing ( |141| ) and (|145|) , we see that the 
leading renormalons dominate the even coefficients of C 77 beginning at the a 4 term, but 
they do not dominate the odd coefficients until much higher order in a s . The reason for 
this is that the coefficient in square brackets in ( |143j ) is 1.28 + (— l) n 1.16. The asymptotic 
coefficients for C 77 are therefore positive and grow like (n — l)!/?^ -1 . However, for odd 
n, there is a near cancellation between the contributions of the renormalons at u = ±| 
that delays the onset of the asymptotic behavior. At preasymptotic values of n, the odd 
coefficients may be dominated by the next higher renormalons at u — ±1. 

We now consider the higher order radiative corrections to the decay rate for J/ip — > 777. 
The large- Nf limit of the short-distance coefficient C 37 (a s ) has not yet been calculated. 
Only the first term in the perturbation series is known ||23|| : 



C 37 (ct s ) = -6.3a + ■ ■ • . (146) 

However we have determined the residues of the u = ±| renormalons, and we can use that 
information to deduce the asymptotic behavior of the coefficients in the large- iVy limit. The 
contributions from the u = — | and u — | renormalons are 



C^'(a s ) —Y 
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14vr 2 -87 a , „,„ 8 e c / 2 m 
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{n-l)\[2foa.{p)\ n . (147) 
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In the MS scheme (C = —5/3) with \i — m, the first few terms in the series are 

C 37 (a s ) 16a - 62a 2 - 660a 3 - 7600a 4 - 160000a 5 - 3100000a 6 + ••• , (148) 

where a = a s (m)/7r. The asymptotic coefficients are almost an order of magnitude large 
than those in ( |144| ) and ( |145| ). The large coefficients are caused by the large residue of the 
u — +2 renormalon. 

One can attribute the large residue of the u = +^ renormalon in C 37 (o; s ) to a poor 
choice of the renormalization scale fi. The contribution from the u — +\ renormalon can be 
decreased at the expense of increasing the contribution from the u = — ~ renormalon. The 
best balance between the two sources of n\ growth is to choose /i so that the two renormalons 
contribute equally to the asymptotic coefficients. This prescription gives the rather small 



value w = 0.28m for the scale. With this choice of scale, the even coefficients in ( |148|) vanish 
and the odd coefficients are decreased in magnitude by about a factor of 2. This scale is 
rather close to the BLM scale for the closely related process J/ip — > ggg, which is fi = 0.31m. 
This suggests that this may indeed be an appropriate choice for /x in the MS scheme. 

We now consider briefly the phenomenological implications of our calculations for char- 
monium decays, where the appropriate coupling constant is a s (m) ~ 0.35. The first terms in 
( |140|) , ( |141| ), and (|146|) give the complete order-a s corrections to the decay amplitudes. They 
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represent radiative corrections of —30% for J/ip — > e + e~ , —19% for r\ c — > 77, and —70% for 
J/ijj — > 777. The order-a^ correction is known only in the large- Nf limit for J/ip — > e + e~ 
and r] c — > 77, and it gives corrections of —7% and +6%. For J/^ — > 777, we only know the 
contribution to the order-a:j; correction from the renormalons at u = ±|, which give a cor- 
rection of —70%. We conclude that the perturbation series for the amplitudes for the decays 
J/ijj — > e + e~ and ?] c — > 77 appear to be well-behaved through next-to-next-to-leading order 
in a s . However, the perturbation series for J/ijj — > 777 appears to be hopelessly divergent. 
One might be able to decrease the order-a;? correction by changing the renormalization scale 
to /i pa m/3, but only at the cost of increasing the order-o; s correction to well over 100%. 



VI. CONCLUSION 

In this paper, we studied the large-order asymptotic behavior of the perturbation series 
for short-distance coefficients in the NRQCD factorization formalism. We calculated the 
Borel transforms of the short- distance coefficients for the electromagnetic decays J r /ip — > 
e + e~ and i] c — > 77 in the large- Nf limit. We found that there exists a universal resummable 
renormalon associated with the cancellation of the Coulomb singularity in the short- distance 
coefficients. We verified that the ambiguities in the short- distance coefficients due to the 
first infrared renormalon are canceled by ambiguities in the long-distance matrix elements 
that contribute through relative order v 2 . We used our calculation to make estimates of 
the higher order radiative corrections to the decay rates for J/ip — > e + e~, r\ c — > 77, and 
J/ip — > 777. This work represents a first step toward studying the effects of renormalons in 
the NRQCD factorization formalism. We hope that further work in this direction can lead 
to an understanding of how to deal with the large radiative corrections encountered in some 
decay rates, such as J/ip — > 777 and J ' /ip — > ggg. 
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FIGURES 

Fig. 1. Tree- level diagram for cc — > e + e~. 
Fig. 2. Tree- level diagrams for cc — > 77. 

Fig. 3. One-loop diagram for the heavy-quark self-energy, with a chain of light quark 
bubbles inserted into the gluon propagator. 

Fig. 4. One-loop vertex correction for cc — * e + e~, with a chain of light quark bubbles 
inserted into the gluon propagator. 

Fig. 5. One-loop diagrams for cc — > 77, with a chain of light quark bubbles inserted into 
the gluon propagator. For each diagram, there is a corresponding diagram with the photon 
lines interchanged. The blob in Figure 5e represents the mass counterterm. 
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